POLYMORPHISM, PHONON DYNAMICS AND CARRIER-PHONON 

COUPLING IN PENTACENE 



Raffaele G. Delia Valle,^ Aldo Brillante,^ Luca Farina,^ Elisabetta Venuti,^ Matteo Masino,^ and Alberto Girlando^'0 

^Dip. Chimica Fisica ed Inorganica and INSTM-UdR Bologna, 
viale Risorgimento 4-y Bologna University, 1-4-0136 Bologna, Italy 
^Dip. Chimica G.I.A.F. and INSTM-UdR Parma, 
CO ^ Parco Area delle Scienze 17/A, Parma University, 1-43100 Parma, Italy 

O . (Dated: February 2, 2008) 

^SJ ' The crystal structure and phonon dynamics of pentacene is computed with the Quasi Harmonic 

Lattice Dynamics (QHLD) method, based on atom-atom potential. We show that two crystalline 
phases of pentacene exist, rather similar in thermodynamic stability and in molecular density. The 
two phases can be easily distinguished by Raman spectroscopy in the 10-100 cm-1 spectral region. 
We have not found any temperature induced phase transition, whereas a sluggish phase change to 
the denser phase is induced by pressure. The bandwidths of the two phases are slightly different. 
The charge carrier coupling to low-frequency phonons is calculated. 
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I. INTRODUCTION 

Despite recent progresses |l| , the development of electronic devices based on organic semiconductors remains an open 
challenge. Considerable attention has been recently devoted to pentacene, due to the claim of mobilities comparable 
to crystalline silicon, and other extraordinary properties 2|. Such claims have been show to be invalidated by scientific 
y misconduct Q , yet pentacene still holds the record of highest room temperature mobility among materials for organic 
Q , thin-film transistors. This record mobility of about 1.5 cm^V~^s~^ has been reached by continuous improvements 
^ in the sample chemical purification 0. On the other hand, the crystal and thin film growth also have paramount 
~~' , importance for device performance, together with a proper understanding of the motion of charge carriers inside the 
^ • crystal. By combining computational and experimental approaches, we have addressed the above problems for the 
specific case of pentacene. Here we review the main results reached so far. 
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II. CRYSTAL STRUCTURES 



The polymorphism of pentacene, both in bulk and in thin film, has been object of intense experimental ^ |H IS 
and theoretical 0, studies. We have analyzed the five complete X-ray structures of pentacene single crystal crystal 
^ reported so far by computing the " inherent structures" , defined as the local minima of the potential energy hypersur- 

d ' face [9]. Each structure was modeled starting from its experimental molecular arrangements and a common ab-initio 
[ molecular geometry. The identity of molecular geometries in all the configurations is a necessary prerequisite for allow- 
I ing identical phases to map into identical potential minima. The theoretical structures of mechanical equilibrium were 
' O determined by minimizing the potential energy, The unit cell axes, angles, molecular orientations and positions 
^ I were freely varied in the minimization. At this stage, we adopted the rigid molecule approximation, retaining the 
molecular geometry fixed. The intermolecular potential was represented by an atom-atom Buckingam model, with 
Williams parameters set IV, plus electrostatic interactions described by the ab-initio charges residing on the atoms 
p. 

With the above described method, we have been able to show that the five reported crystal structures correspond 
5J] to two different "inherent structures", i.e., two local $ minima. One of the minima, labeled as phase C, corresponds 
to the structure determined by Campbell et al. The other polymorph, phase H, corresponds to the structure 
found in the more recent measurements 0. Both polymorphs belong to the triclinic system, space group PI {C}) 
with two molecules per unit cell arranged in the usual herringbone pattern. Fig. ^ compares the two structures, 
which appear to be very similar, the main difference being a shift of the two molecules along the long molecular axis, 
yielding different d(001)-spacings (14.5 vs 14.1 for C and H phases, respectively ^). 
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FIG. 1: Comparison of pentacene C and H polymorphs. 



We have extended the analysis by systematically sampling the potential hypersurface to gain information on the 
global stability of the minima. In the case of pentacene, a quasi-Monte Carlo search yielded several hundred of 
distinct minima, with surprising variety of structural arrangements [Tq| . On the other hand, deep minima are easily 
accessible since their attraction basin tends to be wide. We found that the two deepest minima indeed correspond to 
the experimental C and H polymorphs, which are described correctly. Other deep minima with layered structures, 
which might correspond to the thin film polymorphs found to grow on substrates, were also found. 

After the above analysis of the potential energy <i>, we have accounted for the effects of temperature T and pressure p 
by calculating the structures of minimum Gibbs energy G{p, T) with quasi harmonic lattice dynamics (QHLD) methods 
In this method the vibrational contribution to the Gibbs energy is estimated in the harmonic approximation: 



Here V is the molar volume, and Wq.j are the harmonic phonon frequencies, calculated by diagonalizing the dynamical 
matrix obtained from the second derivatives of with respect to the displacements of the j-th molecular coordinate 
with wavevector q. The second term in Eq. (1) is the zero-point energy, and the last term is the entropic contribution. 
The sums are extended to all phonon frequencies. Since pentacene exhibits ab-initio intra-molecular frequencies as 
low as 38 cm""'^, we had to allow for the coupling between lattice and intra-molecular vibrations, found to be necessary 
in similar cases. We adopted an exciton-like model |^, where the interaction between different molecular coordinates 
is mediated by the intermolecular potential which, being a function of the interatomic distance, depends directly on 
the atomic displacements. Since these correspond to the Cartesian eigenvectors of the normal modes of the isolated 
molecule, we used the ab-initio eigenvectors and the scaled ab-initio frequencies. Intra-molecular modes above 300 
cm~^ were not taken into account, because coupling is important only for low frequency modes. 

Comparison between minimum <i> and minimum G structures shows that the latter, which include vibrational 
effects, brings the calculated volumes from about 4% to about 2% from experimental values, or better. Moreover, 
the minimum G calculations reproduce correctly the thermal expansion of the H phase. The experimental volume 
expansion from 90 to 293 K is 3.5 %, to be compared with a calculated value of 2.7. The coupling between lattice 
and intra-molecular vibrations is important, since the calculated expansion falls down to 1.9 % by neglecting it |^. 

For both phases, the QHLD calculations exhibit a mechanical instability at T > 550 K, where no minimum of 
G{p, T) can be found and the crystal structure diverges. Experimental melting temperatures are given between 540 
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FIG. 2: Typical Raman spectra of pentacene crystals (from Ref. [T^'l. 

and 580 K. As discussed elsewhere loss of mechanical stability in QHLD calculations is not necessarily coincident 
with melting, but is often close to it. A good agreement is also found for the experimental sublimation heat A^ubH 
— 184 ± 10 kJ/mol, which is to be compared with the Gibbs energy calculated at K, Go 174 kJ/mol for both 
phases. 

III. LATTICE PHONON SPECTRA 

Although the above QHLD calculations clearly indicate that both C and H polymorphs are mechanically and 
thermodynamically stable, one might still harbor doubts about the existence of the form reported by Campbell 
et al. 4], since all the recent X-ray studies have failed to reproduce it. Both the minimum $ and minimum G 
calculations indicate that the lowest frequency inter-molecular (lattice) vibrations should differ considerably in C and 
H polymorphs. We have then chosen to use Raman scattering in the lattice phonon region (10-150 cm^^). This 
choice was also supported by the possibility of interfacing optical microscopy to Raman spectroscopy, achieving a 
spatial resolution of about 1 /zm, which enables a careful mapping of the physical features of each crystal sample. 
The spectrometer used was a Jobin-Yvon T64000 triple monochromator, interfaced to an Olympus BX40 microscope, 
and with excitation from a Krypton laser at 752.5 nm. Pentacene crystals (generally metal-like bluish platelets) from 
different commercial sources were used, without further purification, as samples for the Raman measurements. A 
different morphology, i.e., tiny blue micro crystals, was obtained for samples vapor-grown by fast sublimation in Ar 
or N2 atmosphere. The Raman spectra were measured by using the above described different crystal samples. The 
crystals 'physical purity' was checked by mapping their phonon spectra in several spots of the sample. 

Fig. [3 reports typical Raman profiles in the region 20-250 cm^^, where all the lattice, as well as a number of low- 
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frequency intra-molecular modes, are present. Indeed, as we have already pointed up, some intra-molecular modes 
occurring in this spectral region are coupled with lattice phonons. A detailed assignment of the spectra is reported 
elsewhere ^M- Here, we focus on the lowest frequency bands (below 100 cm~^) which differ significantly in 
spectrum I and II. In spectrum I, one identifies four bands almost uniformly spread in the considered frequency range. 
In the case of spectrum II, instead, the three lowest bands are clustered in a narrow frequency range around 55 cm"-'^, 
and one very weak band is found near 100 cm^^. This is precisely the frequency pattern we calculated for polymorph 
H and C, respectively. Thus commercial samples, which give spectrum I, all belong to the H phase, whereas the 
samples obtained by fast sublimation in inert atmosphere give spectrum II and correspond to phase C. 

A further remark should be made on Fig. [5] the spectra originating from commercial samples (polymorph H, 
morphologically platelets) evidence in some cases an additional feature at 50 cm"-'^, as shown by the arrow in 
the spectrum (I + II) of the same figure. This band is clearly due to a physical residual impurity of polymorph C 
(morphologically microcrystals), that can be monitored in regions of the commercial specimens extended only few 
microns, as indicated by the decreasing intensity of the band on increasing the crystal area sampled. The fact that a 
residual presence of polymorph C can be found in different amounts on the measured crystals is certainly related to 
the different degree of physical purity of the used commercial products. 

The Raman indications have been subsequently confirmed by direct X-ray measurements 15J, which indeed yield 
two distinct structures, matching the published measurements for forms C and H. However, the Raman tool is 
more versatile, and allows one to analyze the phase homogeneity, which, as discussed below, might have important 
consequences on the measured mobilities. Although identified only in recent years, the phase H is the most easily 
encountered, irrespectively on the method of preparation, either from sublimation or from solution. One therefore 
expects it is the thermodynamically stable form, even if our calculations indicate a almost equivalent stability for the 
two phases. We have then decided to investigate the possible presence of phase transitions as a function temperature 
and of pressure by using the Raman microprobe technique. At the same time, we could verify that the T, p evolution 
of the phonon frequencies is well reproduced by the computed QHLD frequencies [lolll5l|. 

No phase transition was detected by changing T from 400 to 4.2 K 10]. On the other hand, we have foimd that 
phase C starts to irreversibly transform to the denser H phase by just applying a moderate pressure of about 0.2 Gpa 
|15.J . However, the phase transition goes to completion only when pressures up to at least 5 GPa are applied. Although 
increasing pressure is obviously expected to favor the higher density H phase, the transition mechanism cannot be 
explained on the basis of relative densities only. In fact, thermal annealing at p = 0.6 GPa does not help to speed up 
the structural change started at 0.2 GPa. Together with the sluggish evolution of the transition at room temperature, 
this is an indication that either a high energy barrier must be overcome, or that the driving force becomes weaker 
under pressure. 



IV. CHARGE CARRIER-PHONON COUPLING 



Despite intensive research, we still lack of a proper description of the charge carriers motion in organic molecular 
semiconductors O. The problem is a difficult one, since bandwidths, electron-electron interactions, polarization 
energies, and electron-phonon coupling are all comparable in magnitude. In this paper we shall simply try to estimate 
the relative magnitude of some of the involved parameters. 




FIG. 3: Definition of pentacene largest hopping integrals 
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There is general agreement that the intrinsic mobihties of ohgoacenes exhibit a hopping-to-band transition around 
room temperature |l7|. Here, we consider only the bandlike regime, describing the electronic structure in a tight- 
binding scheme. We have first calculated the intermolecular hopping integrals between pairs of molecules as the 
splitting of HOMO and LUMO levels from the isolated molecule to the pair by using the INDO/S semi-empirical 
method 1 8J . The most important integrals, labeled iai, ta2j tdi, and , are shown in Fig. |3| 

The band dispersions are obtained by: 

e±(k) = {ta, + iaj cos(ka) ± 2 {ta, cos[k(a + b)/2] + t^^ cos[k(a - b))/2]} (2) 

The resulting bandwidths for C and H computed structures are reported in Table|TJ The temperature dependence is 
also given. As it usually happens in this kind of calculations the conduction band (CB) is wider than the valence 
band (VB), contrary to what it would be expected, since in oligoacenes hole mobility is larger than electron mobility. 
We believe this is at least in part due to a computational artifact, since in an Hartree-Fock scheme filled orbital are 
better calculated than virtual orbital energies. So we concentrate on relative variations. The lattice contraction by 
lowering temperature brings about 15 % bandwidth increase. More importantly, a similar difference in bandwidth 
exists between polymorph C and H, wider bands being of course found in the denser polymorph H. Therefore, C 
micro-domains embedded in the commonly grown H phase would likely reduce the band mobility, beside the obvious 
effects of defect scattering at the grain boundaries. 



TABLE I: Calculated VB and CB bandwidths as a function of temperature and structure. 



Structure 




T = K 


T = 150 K 


T = 300 K 


C 


Valence Band (VB) 


548 meV 


511 meV 


470 meV 




Conduction Band (CB) 


588 meV 


548 meV 


508 meV 


H 


Valence Band (VB) 


632 meV 


585 meV 


536 meV 




Conduction Band (CB) 


672 meV 


624 meV 


576 meV 



Extensive theoretical and experimental studies have been devoted to the characterization of charge carrier coupling 
to molecular vibrations, which modulate on-site energies (Holstcin coupling) 19] . On the other hand, very little 
is know about the coupling to inter-molecular, or lattice, phonons (Peierls coupling), despite the fact that lattice 
phonons are thought to affect mobilities even at low temperatures |l6j |. In molecular crystals like pentacene, lattice 
phonons are expected to couple to charge carriers mainly by modulating the amplitude of inter-molecular hopping 
integrals. The corresponding linear Peierls coupling constants are defined by: 

where Qqj is the normal coordinate for the j-th phonon with wavevector q. We have assumed the low-frequency 
phonons, mainly inter-molecular in character, as dispersionless, and computed the Peierls constants for the q = 
eigenvectors only. Coupling to acoustic phonons, which is zero at q = 0, is not considered. Within these approxi- 
mations, symmetry arguments show that only the totally symmetric (Ag) phonons can be coupled to charge carriers. 
To evaluate the Peierls constants we have calculated the hopping integrals, as detailed above, for the equilibrium and 
for geometries displaced along the QHLD eigenvectors. The .g(q,j) are then obtained by numerical differentiation. In 
Table IhI we report the total Peierls coupling constants (sums of the modulations of the four hopping integrals defined 
in Fig. for the H phase at K. 

Tabled shows that several phonons are appreciably coupled to the charge carriers, and that some difference exists 
between the coupling in the VB and in the CB. The overall coupling strength, as measured by the lattice relaxation 
energy, i?LR = J2j9j/'^j^ 35.3 and 26.9 meV for the VB and the CB band, respectively. The largest contribution 
to i?LR (~ 70 %) comes from the two lowest frequencies modes, which are described as librations around the axes 
perpendicular to the long molecular axis, i.e., approximately parallel to the conducting ab plane (Figs. ^andlJJl. 
Given the low frequency of these librations, they are likely to affect mobilities even at low temperatures. A more 
detailed analysis on the pentacene mobilities and their temperature dependence is in progress. 
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TABLE II: Low-energy Ag phonons and Peierls coupling constants of H pentacene at K. 



I 1 

w/cm 


gve/meV 


5CB/meV 


1 -1 
Lo/cm 


ffvB/meV 


gce/meV 


27 


8.2 


7.2 


162 


1.5 


1.6 


61 


7.1 


4.9 


203 


7.2 


5.5 


70 


0.5 


3.1 


206 


2.8 


3.0 


96 


3.7 


5.5 


238 


3.7 


1.5 


132 


1.1 


1.0 


241 


1.0 


1.5 


148 


6.9 


0.0 


259 


5.4 


2.8 


156 


3.8 


5.9 


261 


3.8 


2.4 



Acknowledgments 

Work supported by the Consorzio Interuniversitario Nazionale per la Scienza e Tecnologia dei Material! (I.N.S.T.M.) 
under the Prisma 2002 project, and by the Italian Ministero dell' Istruzione, dell'Universita e della Ricerca (M.I.U.R.). 
We acknowledge enlightening discussions with Karl Syassen. 



[1] C. D. Dimitrakopoulos and D.J. Mascaro, IBM J. Res. & Dev. 45, 11 (2001). 
[2] J. H. Schon, Phys. Stat. Sol. (b) 226, 257 (2001), and refererences therein. 

[3] See the Report of the Investigation Committee on the possibility of scientific misconduct in the work of H. Schon and 

coauthors, available at www.luccnt.com . 
[4] R. B. Campbell, J. M. Robertson, and J. Trotter, Acta Cryst. 15, 289 (1962). 
[5] D. Holmes, S. Kumaraswamy, A. J. Matzger, and K. P. VoUhardt, Chem. Eur. J. 5, 3399 (1999). 

[6] C. C. Mattheus, A. B. Dros, J. Baas, G.T. Oostergetel, A. Meetsma, J. L. de Boer, and T. M. Palstra, Synth. Metals 138, 
475 (2003). 

[7] C. C. Mattheus, G. A. de Wijs, R. A. de Groot, and T. M. Palstra, J. Am. Chem. Soc. 125, 6323 (2003); J. E. Northup, 

M. L. Tiago, and S. G. Louie, Phys. Rev. B 66, 121404(R) (2003). 
[8] E. Venuti, R. G. Della Valle, A. Brillante, M. Masino, and A. Ghlando, J. Am. Chem. Soc. 124, 2128 (2002). 
[9] M. Masino, A. Girlando, R. G. Delia Valle, E. Venuti, L. Farina, and A. Brillante, Mat. Res. Soc. Symp. Proc. 725, 149 

(2002) . 

[10] R. G. Della Valle, E. Venuti, A. Brillante, and A. Girlando, J. Chem. Phys. 118, 807 (2003). 

[11] A. Girlando, M. Masino, G. Viscntini, R. G. Delia Valle, A. Brillante, E. Venuti, Phys. Rev. B 62, 14476 (2000). 

[12] R. G. Delia Valle and E. Venuti Phys. Rev. B 58, 206 (1998). 

[13] A. Brillante, R. G. Della Valle, L. Farina, A. Girlando, M. Masino, and E. Venuti, Chem. Phys. Lett. 357, 32 (2002). 

[14] R. G. Della Valle, E. Venuti, L. Farina, A. Brillante, M. Masino, and A. Girlando, J. Phys. Chem. B, submitted. 

[15] L. Farina, K. Syassen, A. Brillante, R. G. Delia Valle, E. Venuti, and N. Karl, High Press. Res., in press; L. Farina, A. 

Brillante, R. G. Della Valle, E. Venuti, M. Amboagc, and K. Syassen, Chem. Phys. Lett. 375, 490 (2003). 
[16] M. Pope and C. E. Swcnberg, Electronic Processes in Organic Crystals and Polymers, Second Ed., New York, Oxford 

University Press (1999). 

[17] Y. C. Cheng, R. J. Silbey, D. A. da Silva Filho, J. P. Calbert, J. Cornil, and J. L. Bredas, J. Chem. Phys. 118, 3764 

(2003) , and references therein. 

[18] M. C. Zerner, G. H. Loew, R. F. Kirchner, and U. T. Mueller- Westerhoff, J. Am. Chem. Soc. 102, 589 (1980). 
[19] V. Coropceanu, M. Malagoli, D. A. da Silva Filho, N. E. Gruhn, T. G. Bill, and J. L. Bredas, Phys. Rev. Lett. 89, 275503 
(2002), and references therein. 



